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ABSTRACT 


A  computer  program  for  a  transient  one-dimensional  math¬ 
ematical  model  of  combustion  in  porous  medium  was  used  to 
investigate  the  effects  of  several  parameters  on  system 
behavior  .  Results  show  the  effects  of  medium  thickness  , 
permeability  and  reaction  order  on  combustion  .  In  addi¬ 
tion,  the  effect  of  heat  input  on  combustion  was  also 
studied  . 


TABLE  OF  CONTENTS 


I'S- 

'V 


M 


ll 


f 

I 


.J 


»l 


(  ’ 

■'i' 

•  h' 


■i:? 

'iV 

I  .'t 


I A 


INTRODUCTION  .  15 


EFFECT  OF  HEAT  INPUT  ON  COMBUSTION . 19 


A.  INTRODUCTION  AND  DESCRIPTION  .  19 


B .  PROCEDURE 


C .  RESULTS 


1.  CASE  II- 1  SQ  =  18000  BTU/ft.sq  hr  ....  22 

2.  CASE  II-2  SQ  =  20000  BTU/ft.sq  hr  ....  23 

3.  CASE  II- 3  SQ  =  30000  BTU/ft.sq  hr  ....  28 

4.  CASE  II-4  SQ  =  40000  BTU/ft.sq  hr  ....  34 

SUMMARY . 40 


1.  Power  Relation 


2.  Relation  of  Temperature  and  Oxygen 
Concentration  . 


3.  Equilibrium  Temperature  for  Combustion 


EFFECT  OF  THICKNESS  ON  COMBUSTION 


INTRODUCTION 


B .  PROCEDURE 


RESULTS 


1.  CASE  III-l  Thickness  =  0.25  inches 


CASE  III-2  Thickness  =  0.50  inches 


CASE  III -3  Thickness 
CASE  III-4  Thickness 


=  0.75  inches 


=  1.00  inches 


CASE  III-5  Thickness  =  2.00  inches 


CASE  III- 6  Thickness  =  4.00  inches 


CASE  III -7  Thickness  =  6.00  inches 


SUMMARY 


1.  Power  Relation 


Combustion  Start  Time 


IV. 


PERMEABILITY . 89 

A.  INTRODUCTION . 89 

B.  PROCEDURE . 91 

C.  RESULTS . 92 

1.  CASE  IV- 1  Permeability  =  0.00015  ft^ 

(d=D=0.0025  in)  . . 92 

2.  CASE  IV-2  Permeability  =  0.00058  ft^ 

(d=D=0.0050  in)  . . 93 

3.  CASE  IV-3  Permeability  =  0.00233  ft^ 

(d=D=0.0100  in)  . . 97 

4.  CASE  IV-4  Permeability  =  0.00933  ft^ 

(d=D=0.0200  in)  . . 100 

D.  SUMMARY . 103 

1.  Power  Relation . 103 

2.  Combustion  Speed  .  105 

V,  REACTION  ORDER . Ill 

A.  INTRODUCTION . Ill 

B.  PROCEDURE . Ill 

C.  RESULTS . 112 


1.  CASE  V-1  Reaction  Order  =0.25  ....  112 

2.  CASE  V-2  Reaction  Order  =0.50  ....  114 

3.  CASE  V-3  Reaction  Order  =0.75  ....  117 

4.  CASE  V-4  Reaction  Order  =1.00  ....  121 


D.  SUMMARY . 123 

1.  Power  Relation . 123 

VI.  CONCLUSIONS . 132 

LIST  OF  REFERENCES . .  .  134 

INITIAL  DISTRIBUTION  LIST  .  135 


LIST  OF  TABLES 


I 

VARIED  PARAMETERS  DURING 
18000  BTU/FT.SQ  HR  .  . 

THE 
•  • 

PROCESS 

•  •  •  • 

FOR 

SQ  = 

•  •  •  « 

.  27 

II 

VARIABLE  PARAMETERS  DURING  THE  PROCESS  FOR  SQ 
=  20000  BTU/FT.SQ  HR  . .  . 

.  29 

III 

VARIED  PARAMETERS  DURING 
30000  BTU/FT.SQ  HR  .  . 

THE 

PROCESS 

FOR 

SQ  = 

.  35 

IV 

VARIED  PARAMETERS  DURING 
40000  BTU/FT.SQ  HR  .  . 

THE 

PROCESS 
•  •  •  • 

FOR 

SQ  = 

.  41 

V 

VARIED  PARAMETERS  DURING 
THICKNESS  =0.25  INCHES 

THE 

PROCESS 
•  •  »  • 

FOR 

.  55 

VI 

VARIED  PARAMETERS  DURING 
THICKNESS  =0.50  INCHES 

THE 

PROCESS 
*  •  •  • 

FOR 

•  • 

.  59 

VII 

VARIED  PARAMETERS  DURING 
THICKNESS  =0.75  INCHES 

THE 

PROCESS 

•  •  •  • 

FOR 

•  * 

.  64 

VIII 

VARIED  PARAMETERS  DURING 
THICKNESS  =  1.00  INCHES 

THE 

PROCESS 

•  •  •  • 

FOR 
•  » 

.  68 

IX 

VARIED  PARAMETERS  DURING 
THICKNESS  =2.00  INCHES 

THE 
«  * 

PROCESS 

•  •  •  • 

FOR 
•  • 

.  73 

X 

VARIED  PARAMETERS  DURING 
THICKNESS  =4.00  INCHES 

THE 

PROCESS 

•  •  •  • 

FOR 
•  • 

.  77 

XI 

VARIED  PARAMETERS  DURING 
THICKNESS  =  6.00  INCHES 

THE 

PROCESS 

•  •  •  • 

FOR 
•  • 

.  82 

XII 

RESULTS  FOR  ALL  THICKNESS  CASES  .  . 

•  • 

.  88 

XIII 

RESULTS  FOR  PERMEABILITY 

• 

. 

•  • 

104 

XIV 

POWER  RELATION  RESULTS 

• 

•  •  •  • 

• 

108 

XV 

COMBUSTION  AND  EXTINCTION  TEMPERATURE  FOR 
REACTION  ORDER  CASES  . 

129 

LIST  OF  FIGURES 


2.1 

2.2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 

2.10 
2.11 
2.12 

2.13 

2.14 

2.15 

2.16 

2.17 

2.18 
2.19 


Temperature  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ  =  46  seconds  . 

Temperature  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ  =  45  seconds  . 

Oxygen  concentration  vs  X/L  and  time  for  SQ  = 
18DD0  BTU/ft.sq  hr  and  TQ  =  46  seconds  .  .  . 

Oxygen  concentration  vs  X/L  and  time  for  SQ  = 
18000  BTU/ft.sq  hr  and  TQ  =  45  seconds  .  .  . 

Reaction  rate  vs  X/L  and  time  for  SQ  =18000 
BTU/ft.sq  hr  and  TQ  =  46  seconds  . 

Reaction  rate  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ  =  45  seconds  . 

Combustion  and  extinction  carbon  temperature 
for  position  X/L  =  0 . 


Temperature  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  TQ  =  35  seconds  .... 

Temperature  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  TQ  =  34  seconds  .... 


Oxygen  concentration  vs  X/L  and  time  for  SQ  = 
2OO0O  BTU/ft.sq  hr  and  TQ  =  35  seconds  .  .  . 


Oxygen  concentration  vs  X/L  and  time  for  SQ 
20DO0  BTU/ft.sq  hr  and  TQ  =  34  seconds  .  . 


Reaction  rate  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  tQ  =  35  seconds  . 

Reaction  rate  vs  X/L  and  time  for  SQ  =20000 
BTU/ft.sq  hr  and  TQ  =  34  seconds  . 


Extinction  and  combustion  carbon  temperature 
for  position  X/L  =  0 . 


Temperature  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  15  seconds  .... 


Temperature  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  14  seconds  .... 


Oxygen  concentration  vs  X/L  and  time  for 
=30000  BTU/ft.sq  hr  and  TQ  =  15  seconds 


SQ 


Oxygen  concentration  vs  X/L  and  time  for  SQ  = 
3OO0O  BTU/ft.sq  hr, and  TQ  =  14  seconds  .  .  . 


Reaction  rate  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  15  seconds  . 


23 

23 

24 

24 

25 

25 

26 
30 

30 

31 

31 

32 

32 

33 
36 

36 

37 

37 

38 


2.20  Reaction  rate  vs  X/L  and  time  for  SQ  =  30000 

BTU/ft.sq  hr  and  TQ  =  14  seconds  . 

2.21  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =  0 . 

2.22  Temperature  vs  X/L  and  time  for  SQ  =  40000 

BTU/ft.sq  hr  and  TQ  =  9  seconds  . 

2.23  Temperature  vs  X/L  and  time  for  SQ  =  40000 

BTU/ft.sq  hr  and  TQ  =  8  seconds  . 

2.24  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

40000  BTU/ft.sq  hr  and  TQ  =  9  seconds  . 

2.25  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

40000  BTU/ft.sq  hr  and  TQ  =  8  seconds  . 

2.26  Reaction  rate  vs  X/L  and  time  for  SQ  =  40000 

BTU/ft.sq  hr  and  TQ  =  9  seconds  . 

2.27  Reaction  rate  vs  X/L  and  time  for  SQ  =  40000 

BTU/ft.sq  hr  and  TQ  =  8  seconds  . 

2.28  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =1 . .  .  .  . 

2.29  Heat  flux  versus  combustion  temperature  (TQ  ) 

rectangular  plot . ^.  .  .  . 

2.30  Heat  flux  (SQ)  vs  combustion  temperature (TQ„) 

Log  log  plot . ^.  .  .  . 

3.1  Temperature  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1100  F  .  .  .  . 

3.2  Temperature  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1090  F  .  .  .  . 

3.3  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.25  in  Initial  carbon  temperature 

=  1100  F . . . 

3.4  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.25  in  Initial  carbon  temperature 

=  1090  F . T . 

3.5  Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1100  F  .  .  .  . 

3.6  Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1090  F  .  .  .  . 

3.7  Temperature  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  960  F  .  .  .  . 

3.8  Temperature  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  950  F  .  .  .  . 

3.9  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.50  in  Initial  carbon  temperature 

3.10  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.50  in  Initial  carbon  temperature 

=  950  F . . . 


.  58 


3.11 

3.12 

3.13 

3.14 

3.15 

3.16 

3.17 

3.18 

3.19 

3.20 

3.21 

3.22 

3.23 

3.24 

3.25 

3.26 

3.27 

3.28 

3.29 

3.30 


Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  960  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  950  F  ....  58 

Temperature  vs  X/L  and  time  for  thickness  = 

0.75  in  Initial  carbon  temperature  =  900  F  ....  60 

Temperature  vs  X/L  and  time  for  thickness  = 

0.75  in  Initial  carbon  temperature  =  890  F  ....  61 

Oxygen  concentration  vs  X/L  and  time  for 

thickness  =  0.75  in  Initial  carbon  temperature 

=  900  F . . . 61 

Oxygen  concentration  vs  X/L  and  time  for 

thickness  =  0.75  in  Initial  carbon  temperature 

=  890  F . T . 62 


Reaction  rate  vs  X/L  and  time  for  thickness 
0.75  in  Initial  carbon  temperature  =  900  F 

Reaction  rate  vs  X/L  and  time  for  thickness 
0.75  in  Initial  carbon  temperature  =  890  F 

Temperature  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  850  F 

Temperature  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  840  F 


.  62 
.  63 
.  65 
.  65 


Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  1.00  in  Initial  carbon  temperature 
=  850  F . ;  .  .  . 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  1.00  in  Initial  carbon  temperature 
=  840  F . . 


66 


66 


Reaction  rate  vs  X/L  and  time  for  thickness 
1.00  in  Initial  carbon  temperature  =  850  F 


.  67 


Reaction  rate  vs  X/L  and  time  for  thickness 
1.00  in  Initial  carbon  temperature  =  840  F 


.  67 


Temperature  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  760  F 


.  69 


Temperature  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  750  F 


.  70 


Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  2.00  in  Initial  carbon  temperature 
=  760  F . . 


70 


Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  2.00  in  Initial  carbon  temperature 
=  750  F . r  .  .  . 


71 


Reaction  rate  vs  X/L  and  time  for  thickness 
2.00  in  Initial  carbon  temperature  =  760  F 


.  71 


Reaction  rate  vs  X/L  and  time  for  thickness 
2.00  in  Initial  carbon  temperature  =  750  F 


.  72 


3.31 

3.32 

3.33 

3.34 

3.35 

3.36 

3.37 

3.38 

3.39 

3.40 

3.41 

3.42 

3.43 

3.44 

3.45 

3.46 

4.1 

4.2 

4.3 

4.4 


Jemperature  vs  X/L  and  time  for  thickness  = 

.00  in  Initial  carbon  temperature  =  690  F  ....  74 

Temperature  vs  X/L  and  time  for  thickness  = 

4.00  in  Initial  carbon  temperature  =  680  F  ....  74 

Oxygen  concentration  vs  X/L  and  time  for 

thickness  =  4.00  in  Initial  carbon  temperature 

=  690  F . T . 75 

Oxygen  concentration  vs  X/L  and  time  for 

thickness  =  4.00  in  Initial  carbon  temperature 

=  680  F . r . 75 


Reaction  rate  vs  X/L  and  time  for  thickness  = 
4.00  in  Initial  carbon  temperature  =  690  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 
4.00  in  Initial  carbon  temperature  =  680  F 

Temperature  vs  X/L  and  time  for  thickness  = 
6.00  in  Initial  carbon  temperature  =  650  F 

Temperature  vs  X/L  and  time  for  thickness  = 
6.00  in  Initial  carbon  temperature  =  640  F 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  6.00  in  Initial  carbon  temperature 
=  650  F . . 


.  .  76 
.  .  76 
.  .  78 
.  .  79 

.  .  79 


Oxygen  concentration  vs  X/L  and  time  for 

thickness  =  6.00  in  Initial  carbon  temperature 

=  640  F . . . 80 


Reaction  rate  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  650  F  ....  80 

Reaction  rate  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  640  F  ....  81 

Thickness  versus  combustion  temperature 

Retangular  plot . 84 

Thickness  versus  combustion  temperature  Log- log 
plot . . T  .  T  .  .  85 

Thickness  vs  time  for  oxygen  concentration  at 

X/L  =1.  to  reach  zero  Retangular  plot . 86 

Thickness  vs  time  for  oxygen  concentration 

reach  zero  log- log  plot . 87 

Geometry  of  a  typical  cell . 90 

Temperature  vs  X/L  and  time  for  permeability  = 

O.0DO15  ft^  initial  carbon  temperatiure  =  760  F  .  .  93 

Temperature  vs  X/L  and  time  for  permeability  = 

O.0OO15  ft^  initial  carbon  temperatiure  =  750  F  .  .  93 

Oxygen  concentration  vs  X/L  and  time  for 

fermeability  =  0.00015  ft'^  initial  carbon 
emperatiure  =  760  F . 94 


4.5 


Oxygen  comcentration  vs  X/L  and  time  for 
permeability  =  0.00015  ft^  initial  carbo 
temperatiure  =  750  F  . 


on 


94 


10 


4.6  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00015  ft*^  initial  carbon  temperatiure  =  760 

4.7  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00015  ft^  initial  carbon  temperatiure  =  750 

F  • . 

4.8  Temperature  vs  X/L  and  time  for  permeability  = 
O.0OO58  ft^  initial  carbon  temperatiure  =  8S0  F  . 

4.9  Temperature  vs  X/L  and  time  for  permeability  = 

0. 00058  ft^  initial  carbon  temperatiure  =  840  F  . 

4.10  Oxygen  concentration  vs  X/L  and  time  for 

permeability  =  0.00058  ft  initial  carbon 
temperatiure  =  850  F . . . 

4.11  Oxygen  comcentration  vs  X/L  and  time  for 

?ermeability  =  0.00058  ft  initial  carbon 
emperatiure  =  840  F  . 

4.12  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00058  ft^  initial  carbon  temperatiure  =  8S0 

4.13  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00058  ft^  initial  carbon  temperatiure  =  840 

4.14  '  emperatura  vs  X/L  and  time  for  permeability  = 

*  0.00233  ft^  initial  carbon  temperatiure  =  9o0  F  . 

4.15  Temperature  vs  X/L  and  time  for  permeability  = 

0.00233  ft^  initial  carbon  temperatiure  =  9/0  F  . 

4.16  Oxygen  concentration  vs  X/L  and  time  for 
permeability  =  0.00233  ft'^  initial  carbon 
temperatiure  =  980  F  . 

4.17  Oxygen  comcentration  vs  X/L  and  time  for 
permeability  =  0.00233  ft  initial  carbon 
temperatiure  =  970  F  . 

4.18  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00233  ft^  initial  carbon  temperatiure  =  980 

4.19  Reaction  rata  vs  X/L  and  time  for  permeability 
=  0.00233  ft^  initial  carbon  temperatiure  =  970 

4.20  Permeability  versus  combustion  temperature 

retangular  plot  . 

4.2]  Permeability  versus  combustion  temperature  log 

log  plot  . 

4.22  Permeability  versus  f"  Retangular  plot  . 

4.23  Permeability  versus  t*  Log-log  plot  . 

5.1  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.25  Initial  carbon  temperature  =  980  F  .  .  .  . 

5.2  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.25  Initial  carbon  temperature  =  970  F  .  .  .  . 
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5.3 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.25  Initial  carbon 
temperature  =  980  F  . 
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5 . 4 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.25  Initial  carbon 
temperature  =970F  . 
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5.5 

Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.25  Initial  carbon  temperature  =  980  F  . 
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5.6 

Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.25  Initial  carbon  temperature  =  970  F  . 
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5.7 

Temperature  vs  X/L  and  time  for  reaction  order 
=  0.50  Initial  carbon  temperature  =  850  F  .  .  .  . 
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5.8 

Temperature  vs  X/L  and  time  for  reaction  order 
=  0.50  Initial  carbon  temperature  =  840  F  .  .  .  . 
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5.9 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.50  Initial  carbon 
temperature  =  850  F  . 
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5.10 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.50  Initial  carbon 
temperature  =  840  F  . 

119 

5.11 

Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.50  Initial  carbon  temperature  =  850  F  . 
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Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.50  Initial  carbon  temperature  =  840  F  . 
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Temperature  vs  X/L  and  time  for  reaction  order 
=  0.75  Initial  carbon  temperature  =  750  F  .  .  .  . 
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Temperature  vs  X/L  and  time  for  reaction  order 
=  0.75  Initial  carbon  temperature  =  740  F  .  .  .  . 
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5.15 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0,75  Initial  carbon 
temperature  =  750  F  . 
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5.16 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0,75  Initial  carbon 
temperature  =  740  F  . 
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5.17 

Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.75  Initial  carbon  temperature  =  750  F  . 
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5.18 

Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  0.75  Initial  carbon  temperature  =  740  F  . 
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5.19 

Temperature  vs  X/L  and  time  for  reaction  order 
=  1.00  Initial  carbon  temperature  =  670  F  .  .  .  . 
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5.20 

Temperature  vs  X/L  and  time  for  reaction  order 
=  1.00  Initial  carbon  temperature  =  660  F  .  .  .  . 
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5.21 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  1.00  Initial  carbon 
temperature  =  670  F  . 
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5.22 

Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  1.00  Initial  carbon 
temperature  =  660  F  . 

127 

5.23  Reaction  rate  vs  X/L  and  time  for  reaction 
rate=  1.00  Initial  carbon  temperature  =  670  F  . 

5.24  Reaction  rate  vs  X/L  and  time  for  reaction 
order  =  1.00  Initial  carbon  temperature  =  660  F 

5.25  Reaction  order  versus  temperature  Retangular 

plot . 

5.26  Reaction  order  versus  temperature  Log- log  plot 

5.27  Reaction  order  versus  temperature  Semi  log  plot 
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I.  INTRODUCTION 


The  goal  of  this  thesis  is  to  study  the  combustion 
behavior  of  porous  graphite  .  The  mathematical  model  devel¬ 
oped  by  VATIKIOTIS  (1)  is  a  transient,  one-dimensional  model 
of  a  porous  medium  consisting  of  spherical  particles  ,  or 
cylindrical  fibers  of  graphite  . 

The  model  covers  several  aspects  :  mass  transfer,  heat 
transfer,  combustion  ,  airflow  ,  and  temperature  dependency 
of  thermophysical  properties.  Geometric  parameters  of  the 
porous  medium  include  porosity  (p)  and  medium  thickness  (L). 
The  porous  medium  has  interconnected  pores  which  permit  air 
flow  through  the  medium.  This  air  flow  was  modeled  by  the 
continuity  equation  and  Darcy's  Law  .  The  porous  medium  is 
characterized  by  a  number  of  parameters  including  porosity  , 
specific  internal  area  (Z),  tortuosity  (t)  and  permeability 
(m) .  These  parameters  are  defined  as  follows.  Given  a  unit 
square  cell  of  dimension  D,  and  particle  diameter  d,  we  have 


p  =  1  -  ,r  /  4  (d/D)2  (1.1) 

z  =  1/2  ,r  d^/D^  (1.2) 

T  =  1.4  (1.3) 

m  =  Sp^/z^  (1.4) 


Porosity,  p,  is  defined  as  the  ratio  of  void  volume  per 
volume  of  unit  all  .  The  specific  internal  area,  Z,  is  the 


ratio  of  internal  surface  area  to  bulk  volume  .  The 
tortuosity  is  defined  as  the  ratio  of  the  length  of  the 
flovrpath  of  a  fluid  particle  to  the  straight  line  distance. 

The  particle's  size  decreases  as  the  carbon  is  consumed 
and  all  geometric  properties  which  depend  on  particle  diam¬ 
eter  are  functions  of  time  and  position  .  This  model 
considers  that  the  carbon  matrix  remains  rigid  as  the 
particle  diameter  decreases  and  thus  porosity  increases  with 
combustion. 

The  magnitude  of  REYNOLDS  number  defines  the  motion  of 
the  fluid  which  can  be  molecular  ,  viscous  or  inertial.  In 
the  case  of  porous  media  ,  the  flow  is  dominated  by  viscous 
and  inertial  effects.  Darcy's  Law  for  the  range  of  REYNOLD'S 
number  where  viscous  effects  dominate  is  the  equation  that 
governs  the  fluid  flow  in  porous  media.  Two  velocities  which 
distinguish  the  flow  through  the  medium  are  filter  velocity 
Q,  and  pore  velocity  u.  Filter  velocity  is  given  by  : 

Q  =  -  m/;i  (  dP/dx  )  (1.5) 

where  u  is  the  viscosity  of  the  fluid  and  P  is  the  pressure. 

The  pore  velocity  u  is  given  by  : 

U  =  Q/p  (1.6) 

The  hypothesis  of  the  Dupuit-Forcheimer  assumption  is 
that  the  local  pore  velocity  is  greater  than  the  filter 
velocity.  The  actual  velocity  in  a  pore  is  a  function  of 
position  within  the  pore  .  The  Dupuit-Forcheimer  relation 
defines  an  average  velocity  in  a  pore. 

The  model  of  combustion  adopted  was  the  model  of  N.  N. 
Semenov  (2).  The  relation  of  reaction  rate  to  temperature 
and  oxygen  concentration,  and  the  interaction  of  heat  gener¬ 
ation  and  heat  transfer  are  fundamental  to  this  model  .  If 
the  heat  transfer  dominates,  extinction  will  occur.  If  the 


heat  generation  dominates  then  combustion  will  result.  The 
reaction  rate  equation  was  adopted  from  Arrhenius's  Law  for 
a  simple  reaction  where  the  rate  depends  on  the  concentra¬ 
tion  of  the  reactants  and  not  on  the  products.  The  heat 
generated  by  an  exothermal  reaction  is  obtained  by  multi¬ 
plying  the  reaction  rate  by  the  heat  of  combustion  . 

A  study  of  the  curve  of  heat  generation  versus  tempera¬ 
ture  presented  by  Frank-Kamenetskii  distinguishes  two  phases 
during  combustion.  During  the  initial  kinetic  phase  ,  the 
rate  of  reaction  and  the  temperature  are  lowest.  During  This 
phase  ,  there  is  an  excess  supply  of  oxygen.  In  the  kinetic 
regime  the  reation  rate  increases  exponentially  with 
increasing  temperature.  The  kinetic  regime  is  followed  by 
the  diffusion  regime.  This  phase  is  distinguished  by  high 
temperatures  and  reaction  rates  ,  and  the  reaction  is 
limited  by  a  lack  of  oxygen. 

It  is  assumed  that  the  chemical  reaction  for  the  combus¬ 
tion  process  produces  carbon  monoxide  and  carbon  dioxide. 
The  ratio  of  the  mass  rates  of  carbon  monoxide  to  carbon 
dioxide  depends  on  the  temperature.  Increasing  temperature 
results  in  an  increase  of  this  CO  to  CO2  ratio. 

The  model  of  heat  transfer  includes  three  basic  mecha¬ 
nisms  :  convection  ,  conduction  and  radiation  .  The  heat 
transfer  equation  used  by  VATIKIOTIS  (1)  also  includes  a 
heat  generation  term  due  to  combustion.  All  thermophysical 
properties  were  treated  as  temperatures  dependent.  These 
properties  include  viscosity,  conductivity  and  density. 
Energy  balance  equations  were  constructed  for  control 
volumes  of  the  porous  medium  and  the  air  flowing  through  the 
porous  medium  . 

A  third  field  equation  of  the  model  was  obtained  by  a 
mass  balance  of  oxygen  .  This  equation  includes  molecular 
diffusion,  convective  transport  and  a  sink  term  to  account 
for  the  depletion  of  oxygen  due  to  combustion. 
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These  equations  resulted  in  a  system  of  coupled,  nonli¬ 
near  transient  field  equation  .  These  field  equations, 
together  with  boundary  and  initial  conditions  ,  define  the 
problem.  They  were  solved  numerically  by  using  the  Galerkin 
finite  element  method.  Details  of  the  formulation  are 
presented  in  reference  1  . 


II.  EFFECT  OF  HEAT  INPUT  ON  COMBUSTION 


A.  INTRODUCTION  AND  DESCRIPTION 

In  this  section  the  effect  of  heat  input  on  the  combus¬ 
tion  behavior  of  porous  medium  was  investigated.  The  porous 
medium  is  subjected  to  a  heat  source  on  its  surface  Two 
parameters  of  main  importance  are  considered:  The  magnitude 
of  the  heat  flux  of  the  source  to  the  porous  medium  (SQ), 
and  the  time  that  the  heat  flux  remains  active  (TQ). 

To  do  this  analysis,  all  parameters  of  the  system  are 
fixed  except  the  time  period  of  heating  (TQ)  and  the  amount 
of  heat  flux  (SQ) . 

The  purpose  of  this  section  is  to  determine  the  relation 
between  SQ  and  TQ  which  will  result  in  combustion.  This  is 
achieved  as  follows: 

For  a  specific  value  of  heat  flux  SQ,  the  program  is  run 
with  different  values  of  TQ  until  combustion  occurs  at  (TQ)^ 
and  extinction  occurs  at  (TQ)g  ,  where^  is  a  small  time 
period  .  Here  we  took^  seconds  . 

TQj,  =  TQg  +  5  (2.1) 

Four  cases  of  heat  flux  were  studied  in  this  section: 

CASE  II. 1  SQ  =  18000  BTU/ft.sq  hr 
CASE  II. 2  SQ  =  20000  BTU/ft.sq  hr 

CASE  II. 3  SQ  =  30000  BTU/ft.sq  hr 

CASE  II. 4  SQ  =  40000  BTU/ft.sq  hr 

The  fixed  parameters  used  in  all  cases  are: 

Ambient  temperature  =80  F 
Ambient  temperature  =  2117  Ib/ft.sq 
Tortuosty  =  1.400 
Filament  diameter  =  0.0004167  ft 


Thickness  of  matrix  laminate  =  0.0004167  ft 
Thickness  of  porous  medium  =  0.02083  ft 
Gas  constant  for  air  =  53.34  Ibf  ft/lbm  R 
Conductivity  of  filament  =  86  BTU/lbm  H.F. 

Specific  heat  of  filament  =  0.703  Ibm/Cf 
Emissivity  of  filament  =  0.90 

Shape  factor  for  int.HF.XFER  coefficient  =  1.00 

Characteristc  length  of  medium  =  1  inch 

Heat  of  reaction  =  14090  BTU/Cf 

Reaction  order  =  0.50 

Stochiometric  ratio  (fuel/air)  =  0.375 

Reaction  coefficient  =  2065000  Ibm/Cf  H 

Activation  energy  coefficient  =  28840  deg  R 

Pressure  differential  across  thickness  =  50  Ib/ft.sq 

Initial  carbon  and  air  temperature  =80  F 

Initial  uniform  oxygen  concentration  =  0.0172  Ibm/ft.cu 

The  equations  used  in  this  program  permit  three  sets  of 
boundary  conditions  that  are  approximations  to  physical 
situations.  The  boundary  conditions  selected  are  the  best 
approximations  of  a  one  dimensional  model.  The  set  of 
boundary  conditions  used  in  this  case  are  typical  of  thermal 
flow  reactor  without  radiation  from  the  boundary  surfaces  . 
for  the  problems  of  this  investigation,  the  following 
boundary  conditions  were  used  : 

-  No  heat  transfer  from  either  end  of  the  porous  medium, 
dt/dx  =0  at  X/L  =  0  and  X/L  =1  (2.2) 


-  Danckwerts'  boundary  conditions  :  Conditions  for  the 
air  temperature  and  oxygen  concentration  were  assumed.  This 


means  : 


X/L  =  0 


(2.3) 


dTJdyi  =  0  X/L  =  1  (2.4) 

Dg  dit,  /ax  =  u  (^  -  X/L  =  0  (2.5) 

(Kg>K^)  dXg/ax  =  -a  £(Tg^  -  X/L  =  1  (2.6) 

Kgaig/ax  =  p.aCau(Ta  -  Too)  X/L  =  1  (2.7) 

aig/ax  =0  X/L  =  1  (2.8) 

Dga^/ax  =  u(0  -  X/L  =  1  (2.9) 

B .  PROCEDURE 

For  each  case  ,  the  amount  of  heat  flux  (SQ)  was  fixed 
and  several  values  of  shut  off  time  (TQ)  were  chosen  until 
values  of  TQ^  and  TQg  were  found  such  that  : 


C .  RESULTS 


The  results  for  each  case  were  obtained  in  both  numer¬ 
ical  and  graphical  form  .  The  behavior  of  carbon  tempera¬ 
ture,  oxygen  concentration  and  reaction  rate  for  combustion 
and  extinction  are  shown  as  surfaces  versus  time  and  posi¬ 
tion  The  program  also  calculates  the  air  temperature  in  the 
porous  medium.  However  the  air  temperature  surfaces  are  not 
shown  as  they  are  quite  similar  ,  in  most  cases  to  carbon 
temperature . 

There  are  three  surface  graphics  for  each  heat  flux 
case : 

-Carbon  temperature  versus  position  and  time 

-Oxygen  concentration  versus  position  and  time 

-Reaction  rate  versus  position  and  time 

In  edition  ,  there  are  two  curves  for  each  heat  flux 
(SQ)  ,  one  showing  the  variation  of  carbon  temperature 
versus  time  for  combustion  and  the  other  showing  the  varia¬ 
tion  of  carbon  temperature  for  extinction  .  These  is  two 
curves  are  for  the  point  x/L=0.  (i.e.  where  the  air  enters 

the  porous  medium)  . 

1.  CASE  II- 1  SQ  =  18000  BTU/ft.sg  hr 

The  values  for  shut  off  time  bounding  extinction  and 
combustion  are  : 


for  combustion  TQ^,  =  46  seconds 
for  extinction  TQ^  =  47  seconds 

During  the  transient  analyses  ,  system  parameters 
varied  .  for  each  analysis  ,  the  maximum  and  minir.ium  values 
of  some  of  these  parameters  provide  some  insight  into  the 
character  of  the  particular  case  .  These  parameters  and 
their  values  are  given  in  Table  I 

The  graphical  results  are  shown  in  Figures  2.1  to 
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2.  CASE  II- 2  SQ  =  20000  BTU/ft.sg  hr 

For  this  case  the  heat  flux  duration  times  for  the 
conditions  of  extinction  and  combustion  are  : 

for  combustion  TQ^,  =  35  seconds 
for  extinction  TQ^  =  34  seconds 

During  the  transient  history  of  this  case  several 
parameters  varied.  The  maximum  and  minimum  of  these  changing 
variables  are  in  the  Table  II  .  It  is  necessary  to  emphasize 
that  these  extrema  results  do  not  occur  at  the  same  time  or 
same  position  (X/L) .  They  are  maximum  and  minimum  values 
obtained  at  different  times  and  position  during  the  tran¬ 
sient  process.  These  values  show  the  characteristc  behavior 
of  this  case  . 

The  behavior  of  carbon  temperature  ,  oxygen  concen¬ 
tration  and  reaction  rate  are  shown  in  Figures  2.8  to  2.13  . 
These  surfaces  are  an  important  characteristic  of  this  case 
because  they  show  the  evolution  of  the  combustion  process  or 
extinction  process  for  each  position  during  the  transient 
time  . 

3 .  CASE  II-3  SQ  =  30000  BTU/ft.sg  hr 

The  values  of  heating  time  that  result  in  combustion 
and  extinction  for  this  value  of  heat  flux  are  : 


for  combustion  TQ^,  =  15  seconds 
for  extinction  TQg  =  14  seconds 

During  the  transient  analysis  ,  some  variation  of 
values  occured  to  several  parameters.  The  limits  of  varia¬ 
tion  ,  maximum  and  minimum  values  ,  are  shown  in  Table  III  . 
The  process  of  combustion  and  extinction  have  particular 
characteristics  that  determine  the  variation  of  the  values 
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TEUPERPTUPC  SURrRCC  TRC^l  G^R’ZE 


Figure  2.8  Temperature  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  35  seconds 

TEnPCRRTURE  SURFRCi  FROn  GRRr3E 


Figure  2.9  Temperature  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  34  seconds 


OXYGEN  CONC.  SUPEFCE  GRPr3E 


Figure  2.10  Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  35  seconds. 

OXYGEN  CONC.  SUREBCE  GRr=r3E 


Figure  2.11  Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  34  seconds. 
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Figure  2.14  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =  0.. 


of  these  parameters.  The  maximum  and  minimum  values  give  a 
general  idea  of  system  behavior  .  These  limiting  values  do 
not  occur  at  same  position  (X/L)  or  same  time  . 

A  better  idea  of  the  development  of  temperature  , 
oxygen  concentration  and  reaction  rate  can  be  obtained  from 
Figures  2.15  to  2.20  These  surfaces  show  the  values  of  these 
parameters  for  each  point  (X/L)  during  the  transient  time  . 

Figure  2.14  shows  that  at  X/L  =  0.  the  transient 
behavior  for  extinction  and  combustion  are  the  same  until 
TQg  =  14  seconds  ,  the  time  the  heat  flux  is  terminated  . 
After  this  time  the  temperature  for  the  extinction  case 
decreases  to  ambient  temperature  ;  while  for  the  combustion 
case  ,  the  temperature  increases  reaching  a  peak  temperature 
of  3200  F  at  28  seconds  and  thereafter  decreases  to  an  equi¬ 
librium  combustion  temperature  of  about  2500  F  after  about 
50  seconds  (see  Figure  2.14). 


4.  CASE  II-4  SQ  =  40000  BTU/ft.sq  hr 

The  heat  flux  of  this  case  produced  combustion  and 
extinction  of  the  porous  medium  at  the  following  times  : 

for  combustion  TQ^  =  9  seconds 
for  extinction  TQ^  =  8  seconds 

During  the  transient  process  ,  the  output  of  the 
program  showed  parameters  th  at  didn't  change  with  time  and 
others  that  did  .  The  variation  of  these  parameter's  values 
show  the  development  of  the  combustion  and  extinction 
phenomenon  for  this  value  of  SQ  .  The  maximum  and  minimum 
values  of  several  parameters  are  shown  in  Table  IV  .  Again 
these  values  do  not  correspond  to  the  same  time  . 

The  characteristics  of  this  case  ,  carbon  tempera¬ 
ture  ,  oxygen  concentration  and  reaction  rate  ,  as  function 


TEMPERATURE  SURERCE  FROM  GRRPIE 


Figure  2.15  Temperature  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr  and  TQ  =  15  seconds. 


TEMPERATURE  SURPPCE  PROM  GRfiPTE 


Figure  2.16  Temperature  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr  and  TQ  =  14  seconds. 


OXYGEN  CONC.  SURFACE  FROn  GRRF3E 


Figure  2.17  Oxigen  concentration  vs  X/L  and  time 
for  SQ  =30000  BTU/ft.sq  hr  and  TQ  =  15  seconds. 


OXYGEN  CONC.  SURFACE  FROfl  GRRF3E 


Figure  2.18  Oxigen  concentration  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr, and  TQ  =  14  seconds. 


Figure  2.20  Reaction  rate  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr  and  TQ  =  14  seconds. 
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of  time  and  position  (X/L)  are  shovm  as  graphic  surfaces  in 
Figures  2.22  to  2.27  .  These  graphical  results  give  an  idea 
about  the  behavior  of  each  point  of  the  transient  history  . 

For  X/L  =  0  ,  both  the  combustion  and  extinction 

cases  show  similar  behavior  to  this  time  equal  to  about  8 
seconds  ,  when  the  heat  flux  is  terminated  .  After  this 
time,  for  the  extinction  case  ,  the  temperature  decreases  to 
ambient  temperature  .  For  the  combustion  case  ,  the  temper¬ 
ature  at  X/L  =  0  increased  to  about  3300  F  ,  and  thereafter 
decreased  to  around  2500  F  (temperature  of  combustion) 
Both  the  extinction  and  combustion  cases  achieve  new  equi¬ 
librium  temperatures  ,  about  2500  F  for  the  combustion  cases 
and  ambient  temperature  for  the  extinction  cases  ,  after  35 
seconds  (see  Figures  2.28  ). 

D .  SUMMARY 

Here  ,  some  observations  are  made  about  the  effects  of 
heat  flux  input  and  duration  of  heat  flux  on  the  combustion 
process . 

1 .  Power  Relation 

For  each  case  ,  a  pair  of  values  (SQ,  (TQ)^,)  was 
obtained  .  A  plot  of  these  point  is  shown  in  Figures  2.29 
and  2.30  .  It  is  observed  that  the  relation  between  SQ  and 
TQ  of  these  graphics  yields  the  approximate  power  relation 
(equation  2.12)  . 


2 .  Relation  of  Temperature  and  Oxygen  Concentration 

During  the  initial  heating  of  the  porous  medium  , 
the  temperature  increases  while  the  oxygen  is  being 
defleted.  After  the  heat  flux  is  removed  ,  either  extinction 
or  combustion  will  proceed.  If  extinction  occurs  ,  the 
carbon  temperature  decreases  to  ambient  temperature  while 
the  oxygen  is  being  restored  to  ambient  concentreition  . 

If  combustion  occurs  ,then  the  oxygen  inside  the 
medium  goes  to  zero  and  remains  there  .  The  only  oxygen 
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OXYGEN  CONC.  SURPRCE  FROfl  GRPr3E 


Figure  2.24  Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  9  seconds. 


OXYGEN  CONC.  SURFRCE  FROM  GRPF3E 


Figure  2.25  Oxigen  concentration  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  8  seconds. 
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(2.12) 


TQj,  =  3.4  10^°/(SQ)2.09 

present  is  at  the  entrance  at  X/L  =  0  .  The  temperature 

continues  to  increase  to  a  maximum  ,  and  thereafter  the 
temperature  decreases  .  Eventually  the  entire  medium  would 
achieve  a  uniform  combustion  temperature  . 

3 .  Equilibrium  Temperature  for  Combustion 

In  all  combustion  cases  ,  the  behavior  was  similar  . 
During  the  initial  combustion  period  ,  the  temperature  of 
the  medium  at  X/L  =  0  rose  between  2000  F  (CASE  II- 1)  to 
3400F  (CASE  II-4)  .  At  this  time  ,  on  the  other  side  of  the 
medium  at  X/L  =  1.  ,  the  temperature  of  the  medium  was 
approximately  1500  F  .  This  during  thus  early  transient 
stage  there  is  aAT  of  1300  F  (for  CASE  II- 1)  and  1900  F 
(for  CASE  II-4)  across  the  medium  .  As  time  proceeds  the 
temperature  at  X/L  =  0  decreased  ,  while  at  X/L  =1.  ,  the 

temperature  increased  .  That  is  ,  with  progressir'  time  , 
the  medium  moves  towards  a  uniform  temperature  .  In  the 
present  case  ,  this  equilibrium  temperature  of  combustion  is 
around  2400  F  to  2500  F  for  all  cases  . 

The  mechanism  of  equilibrium  of  the  temperature 
during  combustion  is  due  to  the  heat  transfer  by  air  .  As 
the  air  flows  through  the  porous  medium  it  transports  the 
heat  from  the  hot  front  part  of  the  medium  to  the  cooler 
back  part  ,  until  all  points  of  the  medium  have  the  same 
temperature  . 


10000  20000  30000  40000  SOOOO  60000  70000 

5Q  ISURfPCC  HCBT  TLUXJ 


MlNJMUfl  TO  rOR  COMBUST  JON 


Figure  2.29  Heat  flux  versus  combustion  temperature  (TQ^) 

rectangular  plot. 


SQ  VERSUS  TQ 


MINIMUM  Tq  FOR  COMBtmiO 


Figure  2.30  Heat  flux  (SQ)  vs  combustion  temperature  (TQ^, ) 

Log  log  plot. 
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III.  EFFECT  OF  THICKNESS  ON  COMBUSTION 
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INTRODUCTION 


The  porous  medium  has  several  geometric  parameters  that 
can  have  direct  influence  on  the  combustion  process.  The 
effect  of  the  thickness  parameter  on  system  behavior  is 
analysed  in  this  section.  All  other  parameters  are  fixed 
except  the  thickness  of  the  medium. 

The  purpose  of  this  section  is  to  determine  the  relation 
of  the  thickness  on  the  minimum  temperature  to  start  combus¬ 
tion  (i.e.,  the  ignition  temperature). 

For  a  specified  value  of  thickness,  the  program  is  run 
with  different  values  of  initial  temperature.  This  initial 
temperature  is  equal  at  all  points  of  the  medium.  The  lowest 
value  of  initial  temperature  that  results  in  combustion  is 
the  ignition  temperature  Tc.  Extinction  occurs  at  Tc-5  , 
whereS  is  a  small  increment  of  temperature.  The  value  offi 
used  in  this  parametric  analysis  was  10  F  . 

Seven  cases  with  different  thicknesses  were  studied  in 


this  section: 


CASE  III-l  THICKNESS  = 


CASE  III- 2  THICKNESS  = 


CASE  III -3  THICKNESS  = 


0,25  inches 
0,50  inches 
0,75  inches 


CASE  III-4  THICKNESS  =  1.00  inches 


CASE  III -5  THICKNESS  =  2.00  inches 


CASE  III- 6  THICKNESS  =4.00  inches 


CASE  III -7  THICKNESS  =  6.00  inches 


Each  case  is  characterised  by  its  thickness.  The  fixed 
parameters  used  in  all  cases  are: 

-  Ambient  temperature  =  80  F 

-  Ambient  pressure  =  2117  Ib/ft.  sq. 

-  Tortuosity  =  1.400 

-  Filament  diameter  =  0.0004167  ft. 


-  Thickness  of  matrix  laminate  =  0.0004167  ft. 

-  Thickness  of  porous  medium  =  0.02083  ft. 

-  Gas  constant  for  air  =  53.34  Ibf.ft/lbm.R 

-  Conductivity  of  filament  =  86  BTU/lbf  HF 

-  Specific  heat  of  filament  =  0.703  Ibm/cf. 

-  Emissivity  of  filament  =  0.90 

-  Shape  factor  for  int.HF.XFER  coefficient  =  1.00 

-  Heat  of  reaction  =  14090  BTU/CF 

-  Reaction  order  =  0.50 

-  Stoichiometric  ratio  (full/air)  =  0.375 

-  Reaction  rate  coefficient  =  2065000  Ibm/CF  hr 

-  Activation  energy  coefficient  =  28.840  deg  R 

-  Pressure  differential  across  thickness  =  50  Ib/ft'’ 

-  Initial  uniform  oxygen  concentration  =  0.0172  Ibm/ft^ 

B .  PROCEDURE 

For  each  case,  the  characteristic  length  of  the  medium 
(thickness)  was  fixed  and  values  of  the  uniform  initial 
temperature  were  chosen  until  the  value  of  temperature  for 
combustion  (T^)  and  temperature  for  extinction  (Tg)  were 
approximately  equal  (T^  w  Tg)  or  (T^  =  Tg  +  10)  (F)  . 

C.  RESULTS 

For  each  case,  the  results  were  obtained  in  both  numer¬ 
ical  and  graphical  form.  The  behavior  of  carbon  temperature, 
oxygen  concentration  and  reaction  rate  for  combustion  and 
extinction  are  shown  as  surfaces  versus  time  and  position. 

The  numerical  results  from  the  program  show  the  varia¬ 
tion  of  several  parameter  values  during  the  transient  anal¬ 
ysis.  The  air  temperature  inside  the  porous  medium,  the 
griphite  temperature  and  the  oxygen  concentration  are 
presented  at  each  position  (X/L)  at  intervals  of  time.  Other 
system  properties  are  also  given  .  These  properties  charac¬ 
terize  the  process  of  combustion  and  extinction  for 
different  thicknesses.  The  maximum  and  minimum  do  not  occur 
at  the  same  point  or  at  the  same  time. 
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The  graphical  results  are  given  in  the  form  of  surfaces. 
These  surfaces  show  the  behavior  of  carbon  temperature, 
oxygen  concentration  and  reaction  rate  versus  time  and  posi¬ 
tion  (X/L)  for  the  combustion  and  extinction  processes. 

1 .  CASE  III-l  Thickness  =  0.25  inches 

The  values  of  the  initial  uniform  carbon  temperature 
bounding  the  extinction  and  combustion  processes  are: 


for  combustion  Tc  =  1100  F 
for  extinction  Te  =  1090  F 


The  development  with  time  of  the  behavior  of  the 
graphite  temperature,  oxygen  concentration  and  reaction  rate 
during  the  combustion  and  extinction  processes  are  shown  in 
Figures  3.1  to  3.6  The  transient  parametric  analysis  provide 
some  insight  into  the  characterers  of  this  particular  case. 
The  maximum  and  minimum  values  of  these  parameters  are  given 
in  Table  V  . 

2 .  CASE  III-2  Thickness  =  0.50  inches 

The  combustion  and  extinction  process,  in  this  case, 
start  with  the  following  uniform  carbon  temperatures: 

for  combustion  Tc  =  960  F 
for  extinction  Te  =  950  F 

For  this  case,  the  graphical  results  for  graphite 
temperature,  oxygen  concentration  and  reaction  rate  are 
shown  in  Figures  3.7  to  3.12  .  The  characteristics  of  the 
changing  parameters  in  this  case,  are  shown  by  the  maximum 
and  minimum  values  in  Table  VI  . 
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OXYGEN  CONC.  SURFACE  FROM  GRRF3E 


Figure 


Figure 


3.3  Oxygen  concentration  vs  X/L  and  t 
for  thickness  =  0.25  in 
Initial  carbon  temperature  =  1100  F. 

oxygen  CONC.  SURFACE  FROM  GRAFSE 


3.4  Oxygen  concentration  vs  X/L  arid  t 
for  thickness  =  0.25  in 
Initial  carbon  temperature  =  1090  F. 


REACTION  RATE  SURPACE  FROM  GRAF3E 


..  > 

Figure’ 3.5  Reaction  rate  vs  X/L  and  time 
for  thickness  =  0.25  in 
Initial  carbon  temperature  =  1100  F. 


REACTION  RATE  SURFACE  FROH  GRRF3E. 


Figure  3.6  Reaction  rate  vs  X/L  and  time 
for  thickness  =  0.25  in 
Initial  carbon  temperature  =  1090  F. 
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RERCTION  RfiTE  SURPRCE  PROM  GRflP3 


Figure  3.11  Reaction  rate  vs  X/L  and  time 
for  thickness  =  0.50  in 
Initial  carbon  temperature  =  960  F. 

RERCTION  RATE  SURFACE  PROM  GRRP3E 


Figure  3.12  Reaction  rate  vs  X/L  and  time 
for  thickness  =  0.50  in 
Initial  carbon  temperature  =  950  F. 
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.  CASE  III -3  Thickness  =  0.75  inches 


During  this  case  the  minimum  uniform  temperature 
that  starts  the  combustion  process  and  the  maximum  uniform 
temperature  for  extinction  are: 


for  combustion  Tc  =  900  F 
for  extinction  Te  =  890  F 

During  the  transient  analysis,  the  system  parameters 
varied.  For  this  case  ,  the  maximum  and  minimum  values  of 
some  of  these  parameters  provide  some  insight  into  the  char¬ 
acter  of  the  case.  These  parameters  and  their  extreme  values 
are  given  in  Table  VII  . 

Figures  3.13  to  3.18  show  the  graphical  results  for 
carbon  temperature,  oxygen  concentration  and  reaction  rate 
for  this  case. 

TEMPERfiTURE  SURPRCE  PROM  GRRr3E 


Figure  3.13  Temperature  vs  X/L  and  time 

for  thickness  =  0.75  in 
0 

Initial  carbon  temperature  =  900  F. 
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TEMPERATURE  SURFACE  FROM  GRRF3E 


Figure  3.14  Temperature  vs  X/L  and  time 
for  thickness  =  0.75  in 
Initial  carbon  temperature  =  890  F. 


OXYGEN  CONC.  SURFACE  FROM  GRAFSE 


Figure  3.15  Oxygen  concentration  vs  X/L  and  time 
for  thickness  =  0.75  in 
Initial  carbon  temperature  =  900  F. 


REACTION  RATE  SURFACE  FROM  GRRF3E 
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Figure  3  18  Reaction  rate  vs  X/L  and  time 
for  thickness  =  0.75  in 
Ir  il  -arbon  temperature  =  890  F. 

A .  Case  III  ^  Thickness  =  1  00  inches 

The  program  was  run  with  several  values  of  initial 
uniform  graphite  temperature  The  temperatures  bounding 

combustion  and  extinction  are: 


for  combustion  Tc  =  850  F 
for  extinction  Te  =  8*40  F 

The  variable  parameters  of  this  case  attain  maximum 
and  minimum  values  as  given  in  Table  VIII  . 

In  this  case,  the  graphical  results  for  carbon 
temperature,  oxygen  concentration  and  reaction  rate  are 
shown  in  Figures  3.19  to  3.24  . 
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■EN  CCNC.  SUREFICE  FROM  GRRF3 


gure  3.24  Reaction  rate  vs  X 
for  thickness  =  1.00  i 
Initial  carbon  temperature  = 
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5 .  CASE  III-5  Thickness  =  2.00  inches 

The  combustion  and  extinction  process,  in  this  case, 
are  characterized  by  the  following  uniform  initial  graphite 
temperatures : 

for  combustion  Tc  =  760  F 
for  extinction  Te  =  750  F 

The  graphite  temperature  and  other  parameters  varied 
during  these  processes  .  The  variation  of  these  parameters 
characterize  the  system.  The  maximum  and  minimum  values  of 
these  parameters  are  in  Table  IX  . 

The  graphical  results  for  carbon  temperature,  oxygen 
concentration  and  reaction  rate  are  shown  in  Figures  3.25  to 
3.30  . 

Tcf^PERflTURE  SURFflCE  FROM  GRPriE 


Figure  3.25  Temperature  vs  X/L  and  time 
for  thickness  =  2.00  in 
Initial  carbon  temperature  =  760  F. 
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TEMPERATURE  SURFACE  FROM  GRRF3E 


Figure  3.26  Temperature  vs  X/L  and  time 
for  thickness  =  2.00  in 
Initial  carbon  temperature  =  750  F. 

OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  3.27  Oxygen  concentration  vs  X/L  and  time 
, for  thickness  =  2.00  in 
Initial  carbon  temperature  =  760  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRflFSE 


Figure  3.28  Oxygen  concentration  vs  X/L  and  time 
for  thickness  =  2.00  in 
Initial  carbon  temperature  =  750  F. 

REACTION  RATE  SURFACE  FROM  GRAiqg£ 


Figure  3.29  Reaction  rate  vs  X/L  and  time 
for  thickness  =  2.00  in 
Initial  carbon  temperature  =  760  F. 
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REFICTION  RATE  SURFACE  FROM  GRAF3E 


Figure  3.30  Reaction  rate  vs  X/L  and  time 
for  thickness  =  2.00  in 
Initial  carbon  temperature  =  750  F. 

6 .  CASE  III -6  Thickness  =  4.00  inches 

The  minimum  value  of  the  initial  graphite  tempera¬ 
ture  to  start  combustion  of  the  porous  medium  and  the 
maximum  value  of  the  initial  graphite  temperature  that 
results  in  extinction  for  this  case  are: 


for  combustion  Tc  =  690  F 


for  extinction  Te 


680  F 


Other  variable  parameters  of  this  case  are  shown  in 
Table  X  .  The  values  of  these  parameters  are  characteristic 
of  the  combustion  and  extinction  processes  for  this  value  of 
thickness . 

In  this  case,  the  graphical  results  show  the 
behavior  of  carbon  temperature  ,  oxygen  concentration  and 
reaction  fate  in  Figures  3.31  to  3.36 
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CASE  III-7  Thickness  =  6.00  inches 
Below  is  given  the  initial  uniform  carbon  tempera- 
combustion  and  extinction: 


for  combustion  Tc  =  650  seconds 
for  extinction  Te  =  640  seconds 

The  variable  parameters  for  this  case  are  given  in 

The  Figures  3.37  to  3.42  show  the  development  of 
temperature,  oxygen  concentration  and  reaction  rate 
i  position  during  the  combustion  and  extinction 


TEMPERATURE  SURFACE  FROM  GRRF3E 


Figure  3.41  Reaction  rate 
'  for  thickness  =  6. 

Initial  carbon  temperatu 


a 


00  in 

re  =  650  F 


®e:'^ctiqn  rrte  surface  fr:m  grrfze 


Figure  3.42  Reaction  rate  vs  X/L  and  time 
for  thickness  =  6.00  in 
Initial  carbon  temperature  =  640  F. 

D .  SUMMARY 

Here,  some  observations  are  made  about  the  effects  of 
thickness  on  the  combustion  process. 

1 .  Power  Relation 

For  each  case  a  pair  of  values  (thickness,  combus¬ 
tion  temperature)  was  obtained.  Four  plots  of  this  data  are 
shown  in  Figures  3.43  to  3.46  . 

Figures  3.43  and  3.44  show  thickness  versus  combus¬ 
tion  temperature  ,  and  Figures  3.45  and  3.46  show  thickness 
versus  the  time  it  takes  for  the  oxygen  concentration  at  X/L 
=  1.  to  reach. 

It  is  observed  that  the  relation  between  thickness 
and  combustion  temperature  ,  figure  3.43  ha.s  a  form  which 
approximates  an  exponentially  decreasing  function. 

The  log-log  plot  (  Figure  3.44  )  is  close  to  a 

linear  form,  which  yields  the  approximate  power  relation. 

The  temperature  is  in  degrees  fahrenheit  and  thick¬ 
ness  is  in  inches  . 
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-0  16554 


(3.1) 


T  =  850  L 

A  comparison  of  the  computer  program's  results  and 
the  formula's  results  are  given  in  Table  III.  The  maximum 
difference  between  the  analytical  results  and  the  power 
equation  results  is  less  than  3  percent  over  the  range  of 
thickness  between  0.25  and  6.00  inches  . 

This  shows  that  the  formula  is  a  good  approximation 
for  these  cases. 

2 .  Combustion  Start  Time 

During  the  initial  phase  of  the  combustion  process, 
the  oxygen  concentration  at  all  points  of  the  porous  medium 
continuously  decreases  to  zero,  for  all  cases,  the  point  of 
the  porous  medium  that  first  reached  zero  was  X/L  =  1.  The 
time  it  taken  for  that  the  oxygen  concentration  at  X/L  =  1. 
to  reach  zero  is  important  because  it  is  related  to  the  time 
that  combustion  starts. 

Figure  3.45  ,  shows  that  this  time  does  not  increase 
linearly  with  the  thickness.  The  log- log  plot  of  these 
values,  Figure  3.46  ,  shows  that  this  log- log  plot  is 
approximately  a  straight  line. 

The  equation  for  this  line  is: 

T*  =  226.8  (3.2) 

This  equation  can  be  used  to  give  an  estimate  for 
the  time  it  takes  for  combustion  to  start. 


THICKNESS  VERSUS  TEMPERATURE 


«  3  4  ft 

THICKNESS  (INCHES) 


Figure  3.43  Thickness  versus  combustion  temperature 

Retangular  plot. 
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IV.  PERMEABILITY 


A.  INTRODUCTION 

The  permeability  or  hydraulic  conductivity  of  a  porous 
medium  is  the  capacity  of  air  flow  through  the  porous 
medium.  The  air  provides  the  oxygen  for  combustion  .  On  the 
other  hand  ,  air  flow  provides  for  convective  heat  transfer 
within  the  medium  .  Thus  air  flow  provides  for  heat  genera¬ 
tion  from  combustion  ,  as  well  as  heat  transfer  by  convec¬ 
tion  .  The  interaction  of  these  two  mechanisms  determines 
whether  combustion  or  extinction  will  occur  .  If  the  heat 
generation  of  combustion  is  greater  than  the  heat  transfer 
of  convection  ,  then  combustion  will  occur  .  If  the  heat 
transfer  mechanism  dominates  ,  then  extinction  will  occur  . 
The  equation  used  for  permeability  in  the  program  is  : 

m  =  0.2  p3  /  (4.1) 

Where  m  is  the  permeability  ,  p  is  the  porosity  and  z  is  the 
specific  internal  area  .  The  expressions  for  porosity  and 
specific  internal  area  are  : 

p  =  1  -  /  6(d/D)3  (4.2) 


z  =  1/2  ,r  (d/D^)  (4.3) 


z  =  1/2  TT  (d^/D^)  (4.4) 

The  equation  4.3  is  for  spherical  particles  and  equation  4.4 
is  for  cylindrical  particles  .  Equation  4.1  is  the  Kozeny 
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of  po 
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Figure  4.1  Geometry  of  a  typical  cell. 

This  equation  fails  for  highly  porous  fibrous  medias  . 

The  tortuosity  ,t  >  is  the  ratio  of  the  length  of  the 
flow  path  for  a  fluid  particle  to  the  straight  line 
distance.  The  tortuosity  depends  on  the  ratio  d/D  .In  this 
program  tortuosity  was  assumed  to  have  a  constant  value  of 
1.4 

Normally  the  interior  geometry  of  a  porous  medium  can  be 
inregular  and  complex  .  In  such  cases  the  internal  structure 
is  difficult  to  discribe  analytically  .  The  program  used  in 
this  investigation  models  the  porous  medium  as  consisting  of 
cylindrical  fibers  or  spherical  particles  in  a  regular  peri¬ 
odic  geometry  as  shown  in  Figures  4.1  . 

During  the  combustion  process  carbon  is  being  continu¬ 
ously  consumed  and  as  a  result,  the  internal  geometry  is 
changing.  The  medium  becomes  more  porous  as  the  carbon 
decreases.  Therefore,  all  properties  which  depend  upon  fiber 
or  particle  diameter  are  functions  of  time  and  position.  The 


1 


model  used  in  the  program  assumed  that  the  carbon  matrix 
remains  rigid  as  the  particles  diameter  decreases  ,  and  thus 
,  the  porosity  increases  with  combustion  , 

B .  PROCEDURE 

To  study  the  effects  of  permeability  on  combustion  ,  all 
parameters  of  the  program  were  maintained  constant  ,  except 
the  permeability  .  In  order  to  change  the  permeability  and 
keep  the  porosity  constant  ,  it  is  necessary  to  keep  the 
(d/D)  ratio  constant  as  seen  in  the  equations  4.3  and  4.6  . 
This  is  accomplished  by  setting  d/D  equal  toa  .  Then  the 
numerator  of  equation  4.1  remains  constant  ,  however  the 
denominator  changes  because  d  =a  D  and  Z  gives  : 

Z  =  1/2  ir(a  dlB^)  =  constant  /  D  (4.5) 

Z  =  1/2  ir(a^  d^/D^)=  constant  /  D  (4.6) 

Thus  ,  the  dimensions  of  d  and  D  were  changed  propor¬ 
tionally  so  that  the  initial  ratio  (d/D)  remained  constant 
in  all  cases  .  Thus  the  porosity  remained  constant  while  the 
permeability  was  different  for  each  case  . 

The  cases  studied  in  this  section  are  : 


CASE  IV- 1 
CASE  IV- 2 
CASE  IV- 3 
CASE  IV- 4 


PERMEABILITY 

PERMEABILITY 

PERMEABILITY 

PERMEABILITY 


0.00015  ft^  (d=D=0.0025  in) 
0.00058  ft^  (d=D=0.0050  in) 
0.00233  ft^  (d=D=0.0100  in) 
0.00933  ft^  (d=D=0.0200  in) 


For  each  values  of  permeability  the  program  was  run  with 
different  uniform  carbon  temperatures  for  all  points  (X/L)  . 
This  procedure  was  repeated  until  we  obtaind  the  minimum 
temperature  that  results  in  combustion  .  This  temperature 
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is  called  the  combustion  temperature  (T^,)  .  Extinction 

occurs  at  temperatures  below  . 

=  Tg  5  (4.7) 

In  this  section  this  difference  (5  )  was  chosen  equal  at 
10  F  then  : 

Tj,  =  Tg  +  10  (F)  (4.8) 

C.  RESULTS 

For  each  case  ,  the  program  provided  the  values  of  the 
carbon  and  air  temperatures  inside  the  porous  medium  ,  as 
well  as  other  parameters  which  varied  during  the  combustion 
process  . 

1.  CASE  IV- 1  Permeability  =  0.00015  ft^  (d=D=Q.0023  in) 
For  this  value  of  permeability  ,  the  maximum,  carbon 
temperature  that  does  not  result  in  combustion  (extinction 
temperature)  and  the  minimum  temperature  that  results  in 
combustion  (combustion  temperature)  are  : 

for  combustion  =  760  F 

for  extinction  Tg  =  750  F 

The  graphical  results  are  shown  in  Figures  4.2  to 
4.7  .  These  graphical  surfaces  show  the  carboa  temperature  , 
oxygen  concentration  and  reaction  rate  at  each  position  and 
time  .  The  form  of  these  surfaces  shows  the  different 
behavior  of  these  three  parameters  during  the  extinction  and 
combustion  processes. 
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TEMPERFITURC  SURFACE  FROM  GRRF3E 


Figure  4.2  Temperature  vs  X/L  and  time 
for  permeability  =  0.00015  ft^ 
initial  carbon  temperatiure  =  760  F. 


TEMPERATURE  SURFACE  FROM  SRRFSE 


Figure  4.3  Temperature  vs  X/L  and  time 
for  permeability  =  0.00015  ft^ 
initial  carbon  temperatiure  =  750  F. 


OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  4,4  Oxygen  concentration  vs  X/L  and  time 
for  permeability  =  0.00015  ft^ 
initial  carbon  temperatiure  =  760  F. 


OXYGEN  CONC.  SURFACE  FROM  GRPFSE 


Figure  4.5  Oxygen  comcentration  vs  X/L  and  time 
for  permeability  =  0.00015  ft^ 
initial  carbon  temperatiure  =  750  F. 


94 


RERCTION  RATE  SURFRCE  FROM  GRPF3E:  I 


Figure  4,6  Reaction  rate  vs  X/L  and  time 
for  permeability  =  0.00015  ft^ 
initial  carbon  temperatiure  =  760  F. 


RERCTION  RATE  SURFRCE  FROM  GRRF3E 


Reaction  rate  vs  X/L  and  time 

2 


Figure  4.7 

for  permeability  =  0.00015  ft 
initial  carbon  temperatiure  =  750  F. 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANOARDS-1963-A 


2.  CASE  IV-2  Permeability  =  0.00058  ft^  (d=D=0.0050  in) 
Here  the  value  of  permeability  is  0.00058  ft^  which 
results  from  a  pore  diameter  of  0.0050  in  .  For  this  case 
the  combustion  and  extinction  temperatures  are  : 


for  combustion  =  850  0  F 

for  extinction  Tg  =  840  0  F 

From  the  numerical  output  of  the  program  ,  three 
graphical  surface  s  were  created  ;  the  carbon  temperature  , 
oxygen  concentration  and  reaction  rate  as  functions  of  time 
and  position  .  These  surfaces  provide  information  of  the 
development  of  these  three  dependent  variables  during  the 
combustion  and  extinction  process  .  These  surfaces  are  shown 
on  Figures  4.8  to  4.13  . 

TEMPERFITURE  SURFACE  FROM  GRBF3E 


Figure  4.8  Temperature  vs  X/L  and  time 
for  permeability  =  0.00058  ft^ 
initial  carbon  temperatiure  =  850  F. 


TEMPERflTURC  SURFACE  FRQri  GRflF3E 


Figure  4.9  Temperature  vs  X/L  and  time 
for  permeability  =  0,00058  ft^ 
initial  carbon  temperatiure  s  840  F. 


OxrGEN  CONG,  surface  FROM  GRAF3E 


Figure  4.10  Oxygen  concentration  vs  X/L  and  time 
for  permeability  =  0.00058  ft^ 
initial  carbon  temperatiure  =  850  F. 


Figure  4.13  Reaction  rate  vs  X/L  and  time 
for  permeability  =  0.00058  ft^ 
initial  carbon  temperatiure  =  840  F. 

CASE  IV-3  Permeability  =  0.00233  ft^  (d=D=0.0100  in 


For  this  case  the  initial  permeability  is  0.00233 
ft^  .  With  this  value  of  permeability  the  combustion  and 
extinction  temperatures  are  : 


for  combustion 
for  extinction 


Tj,  =  980  O.F 
Tg  =  970  O.F 


The  variation  of  the  dependent  variables  during  the 
combustion  and  extinction  cases  is  shovm  by  a  numerical 
output  of  the  program  .  These  results  are  also  shown  in 
three  graphical  surfaces  ;  carbon  temperature  ,  oxygen 
concentration  and  reaction  rate  on  each  position  (X/L)  and 
time  .  These  graphical  results  are  shown  in  Figures  4.14  to 
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TEMPERATURE  SURFACE  FROM  SRAFSE 


Figure  A. 14  Temperature  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  980  F. 


TEMPERATURE  SURFACE  FRCM 


Figure  4.15  Temperature  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  970  F. 


QXrGEN  CONC.  SURFRCE  'RCM  GRR“3C 


Figure  4.16  Oxygen  concentration  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  980  F. 


O/TGEN  core.  SUPFRCE  GP®r3E 


Figure  4.17  Oxygen  comcentration  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  970  F. 


RErrCTION  RATE  SUR“RCE  FRCM  G?Fir3E 


Figure  4.18  Reaction  rate  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  980  F. 


REPCTION  RATE  SURFACE  FRCr^  GRAFSE; 


Figure  4.19  Reaction  rate  vs  X/L  and  time 
for  permeability  =  0.00233  ft^ 
initial  carbon  temperatiure  =  970  F. 
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4.  CASE  IV-4  Permeability  =  0.00933  ft^  (d=D=0.0200  in) 
The  combustion  and  extinction  temperatures  for  this 

O 

case  ,  with  initial  permeability  equal  to  0.00933  ft^  ,  are: 

for  combustion  =  1180  O.F 

for  extinction  Tg  =  1170  O.F 

for  this  case  ,  it  was  not  possible  to  obtain  the  graphical 
surfaces  due  to  limitations  of  the  graphical  program  . 

D.  SUMMARY 

Here  ,  some  observations  about  the  effects  of  perme¬ 
ability  on  combustion  can  be  understood  by  looking  at  all 
the  results  at  one  time  .  Table  XIII  contains  the  results  of 
all  cases  . 

The  Table  lists  the  temperatures  for  extinction  and 
combustion  ,the  permeability  (m)  ,  the  dimension  of  the  unit 
cell  (d=D)  and  time  t*  .  The  t*  time  is  the  time  at  which 
the  oxygen  concentration  at  a  point  reaches  zero  first  . 
That  point  of  the  porous  medium  that  first  reaches  zero 
oxygen  concentration  in  all  these  cases  is  the  the  position 
point  X/L  =  1  . 

1 .  Power  Relation 

The  pair  of  values  (permeability  ,  combustion 
temperature)  for  each  case  was  plotted  on  a  cartesian  graph. 
An  observation  of  this  curve  shows  that  the  combustion 
temperature  increases  with  increasing  permeability.  Figure 
4.20  shows  that  the  combustion  temperature  increases  mono- 
tonically  and  approaches  an  assindoticptate  as  permeability 
increases  .  The  plot  of  these  same  points  on  a  log- log 
graph  ,  Figures  4.21  ,  shows  that  its  form  is  not  quite  a 
straight  line  . 

It  is  observed  that  the  relation  between  combustion 
temperature  (T^)  and  permeability  (m)  yields  the  approximate 
power  relation  (equation  4.9)  : 


A  comparation  of  the  values  of  the  combustion 
temperature  obtained  by  the  program  and  by  the  power  rela¬ 
tion  is  presented  in  Table  XIV  .  This  Table  shows  that  the 
difference  between  these  results  fail  between  1.2  to  3.2 
percent  .  This  result  shows  that  the  above  power  relation  is 
a  good  approximation  to  estimate  the  temperature  which  will 
start  combustion  for  a  given  permeability  . 

2 .  Combustion  Speed 

In  this  section  ,  the  question  of  how  permeability 
affects  the  time  ^  it  takes  for  combustion  to  begin  is 
addressed  .  We  meassure  the  rapidilly  of  the  combustion 
process  by  the  t*  time  previously  defined  .  Recall  that  t* 
is  the  time  at  which  the  oxygen  first  reaches  zero  at  some 
point  in  the  porous  medium  .  In  all  cases  the  position  X/L  = 
1  is  the  first  point  to  achieve  zero  oxygen  concentration  . 
The  plots  of  permeability  (m)  versus  t*  ,  Figures  4.22  and 
4.23  ,  show  that  t*  decays  exponentially  with  increasing 
permeability  and  approaches  an  assymptotic  value  of  t*  =  0 
as  m  goes  to  infinity  . 

A  very  interesting  observation  is  obtained  from  the 
results  of  Table  XIII  .  It  is  noted  that  the  m^^+^^/m^  (for  i 
=  1,2,3)  remains  pratically  constant  (3.84  ,  4.01  ,  4.00) 
while  the  ratio  takes  on  the  variable  values 
l/(i+l)  (1/1.91  ,  1/2.65  ,  1/3.56)  . 

An  attempt  was  made  to  obtain  an  additional  run  with 

O 

the  permeability  equal  to  0.03733  ft*^  ,  however  the  combus¬ 
tion  activity  was  too  rapid  to  achieve  meaningful  results  . 
Combustion  was  achieved  so  rapidly  that  the  gradient  were 
too  large  for  the  program  to  handle  . 

The  permeability  is  a  parameter  that  has  great 
influence  on  the  combustion  process  for  porous  medium  . 


PERMEABILITY  VERSUS  TEMPERATURE 
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Figure  4.20 


Permeability  versus  combustion  temperature 
retangular  plot. 
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V.  REACTION  ORDER 

A.  INTRODUCTION 

The  expression  for  reaction  rate  for  carbon  reacting  in 
air  used  in  the  program  was  the  expression  resulting  from 
the  experimental  data  of  Parker  and  Hotel  . 


Reaction  order  is  the  exponent  of  the  partial  pressure 
of  oxygen  in  the  expression  of  the  reaction  rate  .  The 
program  permits  the  user  to  select  any  value  for  this  param¬ 
eter  .  In  all  of  the  previous  chapters  ,  a  value  of  n  equal 
to  0.5  was  used  .  In  this  chapter  we  investigate  the  effect 
of  n  on  the  combustion  process  . 


B .  PROCEDURE 

In  this  section  ,  all  parameters  are  fixed  to  their 
basic  values  and  only  the  reaction  order  is  changed  in  each 
case  . 

The  program  was  run  ,  for  each  value  of  reaction  order  , 
with  different  value  of  the  uniform  initial  carbon  tempera¬ 
ture  until  we  obtained  the  minimum  temperature  which  results 
in  combustion  .  This  temperature  is  the  combustion  tempera¬ 
ture  (Tj,)  .  The  maximum  temperature  that  results  in  extinc¬ 
tion  is  the  extinction  temperature  (Tg).  they  differ  by  an 
inf initesimalfi  (T^  =  Tg  +5)  •  The  magnitude  offi  was  chosen 
equal  to  10  F  for  this  study  . 

In  order  to  study  the  influence  of  the  reaction  order  on 
the  combustion  process  ,  four  different  values  of  this 
parameter  were  chosen  .  These  cases  are  : 
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CASE  V-1 
CASE  V-2 
CASE  V-3 
CASE  V-4 


Reaction  Order  =  0.25 
Reaction  Order  =  0.50 
Reaction  Order  =  0.75 
Reaction  Order  =  1.00 


The  fixed  parameters  used  in  all  cases  are  : 

Ambient  temperature  =  80  F 
Ambient  pressure  =  2117  Ib/ft.sq 
Tortuosty  =  1.400 
Filamento  diameter  =  0.0004167  ft 
Thickness  of  matrix  laminette  =  0.0004167  ft 
Thickness  of  porous  medium  =  0.02083  ft 
Gas  constant  for  air  =  53.34  Ibf  ft/lbm  R 
Conductivity  of  filament  =  86  BTU/lbm  H.F. 

Specific  heat  of  filament  =  0.703  Ibm/Cf 

Emissivity  of  filament  =  0.90 

Shape  factor  for  internal  HFxFER  =1.00 

Characteristc  length  of  medium  =  1  inchess 

Heat  of  reaction  =  14090  BTU/Cf 

Stochiometric  ratio  (fuel/air)  =  0.375 

Reaction  coefficient  =  2065000  Ibm/Cf  H 

Activetion  energy  coefficient  =  28840  deg  R 

Pressure  differential  across  thickness  =  50  Ib/ft.sq 

Initial  uniform  oxygen  concentration  =  0.172  Ibm/ft^ 


C .  RESULTS 

For  each  cases  the  results  were  obtained  in  both  numer¬ 
ical  and  graphical  form  .  These  results  show  the  variation 
of  some  parameters  during  the  combustion  and  extinction 
processes  .  The  results  obtained  for  each  value  of  reaction 
order  are  presented  below 

1 .  CASE  V-1  Reaction  Order  =  0.25 

With  this  value  of  reaction  order  ,  the  combustion 
and  extinction  processes  define  the  following  characteristic 
temperatures  : 
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for  combustion  T^.  =  980  o.F 
for  extinction  Tg  =  970  o.F 


The  analysis  of  the  numerical  output  of  the  program 
shows  that  some  parameters  varied  during  the  combustion  and 
extinction  processes  .  The  limits  of  these  parameters  show 
the  character  of  the  particular  kind  of  process  which 
occured  . 

The  graphical  results  show  three  surfaces  that 
represent  the  evolution  of  the  carbon  temperature  ,  oxygen 
concentration  and  reaction  rate  at  each  position  during  the 
transient  history  .  These  graphical  surfaces  are  shown  in 
Figures  5.1  to  5.6  . 


TJIMPlRRTURE  SUPFRCi  TROr*  •3RPf3E: 


Figure  5.1  Temperature  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  980  F. 
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TEMPERRTURE  SURFflCE  rROM  GRfVTC 
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Figure  5.2  Temperature  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  970  F. 


OXfGEN  CONC.  SURFACE  GRAFSE 
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Figure  5.3  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  980  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRflF3E 


Figure  5.4  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  970  F. 

PERCTION  RATE  SURFACE  FPO*^  GRAFSE 


Figure  5.5  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  980  F. 


REACTION  RATE  SURFACE  FROM  GRRF3E 


Figuie  5.6  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.25 
Initial  carbon  temperature  =  970  F. 

2 .  CASE  V-2  Reaction  Order  =  0.50 

For  reaction  order  equal  to  0.50  ,  the  initial 

uniform  carbon  temperatures  that  separate  the  combustion  and 
extinction  processes  are  : 


for  combustion  =  980  F 
for  extinction  Tg  =  970  F 

The  program's  numerical  output  included  the  values 
of  some  parameters  that  varied  during  the  transient  process. 
The  variation  of  the  values  characterize  the  combustion  or 
extinction  process.  The  graphical  results  are  shovm  in 
three  surfaces:  carbon  temperature,  oxygen  concentration  and 
reaction  rate  versus  position  (X/L)  and  time.  These  graph¬ 
ical  surfaces  are  shovm  in  figures  5.7  to  5.12  . 


TEMPERflTURC  SURFACE  FROM  GRfiF3E 


Figure  5.7  Temperature  vs  X/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  850  F. 


temperature  Surface  from  grafse 


Figure  5.8  Temperature  vs  X/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  840  F. 


OXrGEN  CONC.  SURrnCE  prom  GRRr3E 


Figure  5.9  Oxygen  concentration  vs  li/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  850  F. 


OXYGEN  CONC.  SURPRCE  PROM  GRRP3E 


Figure  5.10  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  840  F. 


Figure  5.11  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  850  F. 


REACTION  RATE  SURFACE  FROM  GRRF3E 


Figure  5.12  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.50 
Initial  carbon  temperature  =  840  F. 
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TEMPERATURE  SURPPCE  FROM  GRPF3E 


Figure  5.14  Temperature  vs  X/L  and  time 
for  reaction  order  =  0.75 
Initial  carbon  temperature  =  740  F. 

oxygen  CONC.  SURFACE  FROM  GRPF3E 


Figure  5.15  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  0.75 
Initial  carbon  temperature  =  750  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  5.16  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  0.75 
Initial  carbon  temperature  =  740  F. 


REACTION  RATE  SURFACE  FRO^  GRAF3E 


Figure  5.17  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.75 
Initial  carbon  temperature  =  750  F. 


REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure  5.18  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  0.75 
Initial  carbon  temperature  =  740  F. 

4.  CASE  V-4  Reaction  Order  =  1.00 

For  this  case,  the  temperature  of  the  combustion  and 
extinction  processes  are: 


for  combustion  =  670  F 
for  extinction  Tg  =  660  F 


The  Figures  5.19  to  5.24  show  the  behavior  of  the 
carbon  temperature,  oxygen  concentration  and  reaction  at 
each  position  and  time. 

D .  SUMMARY 

1.  Power  Relation 

For  each  case,  a  pair  of  values  (reaction  order, 
combustion  temperature)  was  obtained. The  points  were  plotted 
in  cartesian  coordinates  in  Figure  5.25  .  A  log- log  plot  is 
shown  in  Figure  5.26  ,  and  a  semi  log  plot  is  shown  in 
Figure  5.27  .  It  is  observed  that  the  relation  between 


TD1PDWTURE  SURERCE  PROM  GRRPIE 


Figure  5.19  Temperature  vs  X/L  and  time 
for  reaction  order  =  1.00 
Initial  carbon  temperature  =  670  F. 

reaction  order  (n)  and  combustion  temperature  (T^,)  of  figure 
5.27  yields  the  approximate  power  relation  (equation  5.2)  . 

Table  5.25  shows  a  resume  of  each  case,  combustion 
temperature,  extinction  temperature  and  the  value  of  the 
combustion  temperature  obtained  from  Equation  5.2  and  the 
percentage  error  .  An  observation  of  these  results  leads  to 
the  conclusion  that  this  power  relation  is  a  good  estimation 
for  combustion  temperature. 
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TEMPERATURE  SURFACE  FROM  GRAF3E 


Figure  5.20  Temperature  vs  X/L  and  time 
for  reaction  order  =  1.00 
Initial  carbon  temperature  =  660  F. 

OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  5.21  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  1.00 
Initial  carbon  temperature  =  670  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRflr3E: 


Figure  5.22  Oxygen  concentration  vs  X/L  and  time 
for  reaction  order  =  1.00 
Initial  carbon  temperature  =  660  F. 

REACTION  RATE  SURFACE  FROM  GRAF3E  ^ 


Figure  5.23  Reaction  rate  vs  X/L  and  time 
for  reaction  rate=  1.00 
Initial  carbon  temperature  =  670  F. 


REACTION  RATE  SURFACE  FROM  GRPF3E 


Figure  5.24  Reaction  rate  vs  X/L  and  time 
for  reaction  order  =  1.00 
Initial  carbon  temperature  =  660  F. 


TABLE  XV 


COMBUSTION  AND  EXTINCTION  TEMPERATURE 
FOR  REACTION  ORDER  CASES 

temperatures  (qF) 


extinction 

combustion 

%  error 

reaction 

program 

program 

power 

order 

output 

output 

relation 

0.25 

970 

980 

980 

0.0 

0.50 

840 

850 

863 

1.6 

0.75 

740 

750 

760 

1.4 

1.00 

670 

670 

670 

0.0 
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COMBUSTION  TEMPERATURE  (F) 


REACTION  ORDER  VERSUS  TEMPERATURE 


Figure  5.25  Reaction  order  versus  temperature 

Retangular  plot. 


REACTION  ORDER  VERSUS  TEMPERATURE 


Figure  5.26  Reaction  order  versus  temperature 

Log- log  plot. 
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REACTION  ORDER  VERSUS  TEMPERATURE 


VI.  CONCLUSIONS 


In  each  of  the  previous  chapters,  the  effects  of  a 
single  parameter  on  the  combustion  behavior  of  a  porous 
medium  was  investigated.  In  Chapter  2,  we  considered  the 
effect  of  heat  flux  (SQ)  and  its  duration  (TQ)  on 
combustion.  In  Chapter  3,  the  effect  of  medium  thickness  (L) 
on  behavior  was  investigated.  Chapter  4  was  a  study  of  the 
effect  of  permeability  (m)  on  system  behavior;  and  finally 
in  Chapter  5,  the  effect  of  reaction  order  (n)  was 
investigated . 

The  major  conclusions  from  these  investigations  are: 

-  The  relation  between  SQ  and  TQ  is  not  linear  ,  a  power 
relation  gives  rasonable  approximation. 

-  During  the  combustion  process,  for  the  SQ  cases, 
initially  there  is  a  large  temperature  gradient  across 
the  medium  with  T  (  X/L  =  0.)  >>  T  (  X/L  =1.)  .  With 
passing  time,  this  gradient  decreases  until  there  is  a 
uniform  combustion  temperature  over  the  entire  medium. 
When  this  occurs ,  there  is  no  oxygen  within  the 
medium,  and  the  medium  is  continuosly  consumed  by 
combustion. 

-  The  relation  between  thickness  and  combustion  tempera 
ture  is  not  linear  .  Increasing  the  thickness  results 
in  decreasing  of  combustion  temperature.  If  the 
thickness  increased  with  a  constant  rate  c  (  L^^.^/1^  = 
c  )  the  time  that  marks  the  start  of  combustion  is  not 
constant  ,  it  increases  with  increasing  thickness  ( 
i.e.  t  i+l/t  i  =  b(L)  )  where  b  depends  upon  the 
thickness  . 

-  The  relation  between  permeability  and  minimum  combus 
tion  temperature  is  not  linear. 

-  An  increase  of  the  permeability  shows  that  the  combus 
tion  temperature  increases  to  an  asymptotic  limit. 


-  Permeability  has  a  significant  influence  on  the  speed 
of  the  combustion  process.  When  the  permeability  is 
increased,  the  time  that  marks  the  start  of  the 
combustion  process  decreases.  For  large  values  of 
permeability,  the  combustion  may  occur  in  a  fraction 
of  a  second. 

-  If  the  permeability  is  increased  with  a  constant  rate 

c(  =  c  m^  )  ,  the  time  that  measures  the  start  of 

combustion  decreases  . 

-  The  relation  between  reaction  order  (n)  and  combustion 
temperature  (Tc)  is  not  linear. 

-  An  observation  of  the  graphical  surface;;,  for  all 

cases  of  this  investigation  ,  shows  a  similar  form  for 
all  combustion  and  extinction  cases.  For  all 
combustion  cases,  it  was  noted  that  the  oxygen  within 
the  porous  medium  being  consumed  with  passing  time. 
The  first  point  at  which  the  oxygen  is  depleted  is 
always  X/L  =  1.  With  passing  time  the  point  with 

smaller  X/L  reach  zero  oxygen,  until  all  points  have 
zero  oxygen.  At  this  time  the  only  oxygen  is  at  the 
interface  of  the  porous  medium  and  the  ambient 
environment.  At  this  time  the  combustion  process 
becomes  a  surface  recession  phenomena.  At  the  same 
time  the  oxygen  is  being  depleted  the  temperature  is 
increasing.  Eventually  when  all  the  interior  oxygen 
has  been  consumed,  there  is  a  temperature  gradient 
across  the  medium  with  T^(X/L=0)  >  Tj,(X/L=l). 

This  thesis  concludes  with  two  recomendations  .  First  , 
it  is  recommended  that  additional  parameter  studies  he 
conducted  .  For  example  ,  it  would  be  useful  to  determine 
the  effects  of  pressure  differential  ,  conductivity  ,  and 
other  thermophysical  properties  on  combustion  .  Another 
recommendation  is  to  attempt  to  obtain  a  relation  between 
combustion  temperature  and  several  parameters  . 
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